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INTRODUCTION 



In an academic sense alone, a mathematical investigation 
of the stability of towed bodies and systems is of interest 
only if it adds to the existing knowledge in the general field 
of the course-keeping qualities of vessels. In a more practi- 
cal sense, it is an important and, in fact, necessary primary 
design aid which yields specific criteria to be satisfied for 
stable operation. These criteria are composed of terns which, 
for the most part, have distinct physical significance and 
which can be varied in design studies. Through an appropriate 
analysis of such terms, desirable operating characteristics may 
be found to a considerable degree of accuracy while such a de- 
sign is still in the "drafting board" stage. 

This paper deals with the analysis of a system composed 
of a towboat or tug and a single barge. Both bodies are con- 
sidered to be unrestrained, hence the stability of the system 
will be dependent on the characteristics of the component 
parts. A simple automatic control is specified for the towing 
body in order that the system should have directional rather 
than dynamic stability (or instability). 

In general the methods followed are those of Euler, Routh 
and Hurwitz and deal with systems of simple linear differential 
equations only. Laplace Transformations were not applied, 
since the author has found that many engineers are unfaniliar 
with the properties of this concept . 



MATHEMATICAL DEVELOPMENT 



The subsequent analysis presupposes the following basic 
simplifying assumptions: 

1. Both bodies have freedom in the horizontal plane 
alone, i.e., roll and pitch are ignored. 

2. The towline is assumed horizontal and taut at all 

times. 

3. Angular displacements are considered to be small, so 
that second order and higher terms may be neglected in the 
final equations of motion. 

With these premises in mind, consider a set of right 
hand orthogonal axes fixed relative to the earth and denoted 
by (* 0 .y 0 ). Consider also similar systems of axes- ( z ,y) and 
(x,y) fixed in the barge and towboat respectively. Such a 
system is shown in Fig. 1 together with other symbols neces- 
sary to define completely the orientation of the system at any 
time, A complete and detailed nomenclature is given in Appen- 
dix I. 




Kinematic Relations 



Certain expressions are obtained from the towing geometry. 
These are presented without comment based on fig. 1 of Ref. 1. 
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T x = T Cosf . T y = - T Sinf 

T % = - T CosTj , Ty = T Sin 7 ? 

Differential Equations of Motion 

Since axes fixed in the bodies are essential to the de- 
velopment it is necessary to follow Euler, Ref. 2 in determining 
appropriate mass-acceleration terms. For the purpose of con- 
tinuity these expressions will now be derived. 

Choose a set of axes fixed in space and a set of axes 
fixed in the body as shown in Fig. 1. If P x and Py are the mo- 
mentum expressions relative to the axes fixed in the body then 
(according to Newton) the force components along the axes fixed 
in space are equal to the components of the time rate of change 
of momentum along these axes. Thus 
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where 



l\ = St, u l> y = W s v 

Then 

V = SI, ii Wj i> 0 

V = A/ ? v f l/j a 0 
Similarly 

\ = »/ ( u W ? V 0 

F = ,W 2 V + St, U 0 

The moment of inertia terms are 
/V = / 0 and V = / 0 

We consider here Ref, 1 and Ref. 2 and note that 
,W, = V (1 * k t ) , = St (1 + k t ) 



♦ ; 
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(9) 
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11 ) 



where M is the absolute mass of the ship; l B is the moment of 
inertia of the body about a vertical axis taken through its 
center of gravity, I w is the moment of inertia of the fluid 
displaced by the body taken about a vertical axis through its 
center of gravity and fe,, fc 2 anc * ^ 31,6 appropriate coeffi- 



cients of acce33ion to inertia. Similar expressions may be 
written for the "bar" terms, thu3 

W, = V (1 * k t ) W ? - SI (I ♦ k ) 

(1; 

1 ~ k ' 

We are now in a position to write the differential equa- 
tions of motion: 

1. Barge 



a) 


Longitudinal Force 






SI, u - >/ ? v 0 = A (u,v,0J - T x 


(13) 


b) 


Lateral Force 






W, v + ,l/ 1 u \L = Y (u.v.^i) + t 


rih) 


c) 


Yawing Moment 






/ = V ('ll t>,0) + /Vfj.) 
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Towboat 




a) 


Longitudinal Force 






l/j n v w = A (U.v.ib.b) + T x 


(16) 


b) 


Lateral Force 






,1/, ii + u 0 = jp + Ty 


(17) 


c) 


Yawing Moment 


18) 




7 0= A (u.v.tji.h) * Nf T) 



3. The towboat rudder equations for rudder angle pro- 
portional to heading change are: 

a) Rudder order 

& = « 0 v 19) 

b) Constant time lag 

V = h ♦ T & - ) 

dt 



Linear Simpllfica* ions 



Unle the previous . y' toms are linearized, a cla ical 
solutior io extremely difficult to obtain, in fact in a non- 
linear form the problem can best be treated through the use of 
automatic computing device--. However, only a closed form solu- 
tion will be considered herein, hence the systems will now be 
linearized . 

It is first assumed that V £ V , i.e. that both the tow- 
boat and barge are traveling at the same speed and that both P 
and F take on only small values. Then the following relations 
are obtained from Figure 1 

u = V Cos/ 3 .C V 

V = V Si n P % VP 

u - V Cos 1 £ V 

v = -V Si nP ~ -VP 



Noting these relations, the truth of the following iden- 



tities is obvious 

A (u, v,i I') = A (P U<) 

Y ( u , v ,V>) = Y (P,'l>) 

iV (u.v.V) = N (P,'<) 

Y (Tt.v J.b) = V (P.p.b) 

Y (u ,v,yJi,b) = Y 

N (u.v.h.T) = N (P,J<,~b) 



( 22 ) 



Limiting all displacements to the first order of small 
quant itie°, relat ions (22) may be expanded in Taylor's Series 
to obtain + he following: 
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where the subscript e denotes the equilibrium values about which 
the expansion takes place. Since it is assumed that initially 
the system is on straight course and in line it is obvious that: 



U = < V e = y e = N e = 0 



and 



8 - \!) = 8 =0 =8 - 0 
e e e e e 



In addition, for reasons of symmetry. Ref. 3. 

3A' _ 3i _ 3A _ BA _ 3.Y _ ^ 

3/3 3£ aJr 3^ 38 

Relations (22) then reduce to- 

v = V e Y = YpP + Yyi> 

N = ,Y/ - /V^ 

v = X e 7 = YpP + >>7 + FjS 

V = Sgp - i- VjS 



(2u) 



where the subscript notation has replaced the conventional par- 



tial derivative form. 



Continuing th » linearization to include expre 'iorr 
(2). and (3) the following relation- are established 
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In a similar fashion expressions (d), (9) and (10) become 



\ = 0 
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V = / I/' 

V = 0 

7 = 17, v ~ - u 2 v 0 
v = 7 ~ 



(26) 



With these linearizations the equations of motion become: 
1. Barge 

a) Longitudinal Force 

0 = \ e + 7 (27) 

b) Lateral Force 

V, V' " « 2 »’ ^ = >/ + 
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(26) 



c) Yawing Moment 
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( 29 ) 



2. Towboat 



a) Longitudinal Force 



0 - -T 



(30) 



b) Lateral Force 

(/j V j> - W 2 F t = Ty5 ♦ Fjj ♦ fjo ♦ T [“ 'JyV^-/ 5 ) 



c) Yawing Moment 
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r F + if ( ' p ~ /3)d ' ~ jf ^-^ dt 

(32) 



d) Rudder Order 
I* = a D 0 

e) Constant Time Lag 

r = i + ? «*! 

dt 



(33) 



(3M 



Equations (27) and (JO) represent steady-state conditions 
and need not be considered further in the stability analysis. 
Non-Dimensionalization of Equations of Motion 

Since the hydrodynamic forces may best be expressed in 
terms of dimensionless coefficients, one divides equation 

; U?> b v£/ 3 r 2 



(28) by |f 2 f 2 

(3D by £pV 2 



: and (32) by ££ V 2 



Furthermore the following t ransf ormation is introduced for the 
independent variable * 






Here s i„- f he nurb°r of tow lengths traveled, t the true 

time, L the tow length and l *he linear velocity of the system. 

Since V is con idered substan + ially constant the following 

relation is written — 
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Towboat 

a) Lateral Force 
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STABILITY ANALYSIS 



The equations of notion of the bodies under consideration 
together with the coupling effects due to the towline con- 
straint and the inposition of the control functions on the tow- 
ing body have been presented in the previous section. 

These equations are of general significance within the 
limits of linear theory with the exception that a constant tine 
lag has been assumed for the control function. This assumption 
has been made purely to show that time lags of a sort can be 
included in the development. However in the subsequent analy- 
sis the control time lag s" will be assumed to be zero. No 
loss in generality will occur since it is evident that the in- 
clusion of such a term will only serve to increase the order 
of the stability equation by one. Its exclusion will not make 
the analysis any less valid but will to some degree serve to 
simplify it. 

To remove the integral signs in equations (37), (38). 

(39) and (hO ) . we operate on them by and hence raise 
their order by one. The equations then become: 

1. Barge 



a) Lateral Force 




(u3) 




b) Yawing Moment 
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Towboat 

a) Lateral Force 
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c) Rudder 



S = a Q 0 



(h7) 



These are linear homogeneous equations in the dependent 
variables , Z5 , 0 . , 8, hence their solutions, excluding 

the highly improbable cases where the roots are equal, are 
given by the following equations: 



* = 2 </» e CTls 
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a i s 



It is our aim to determine: 



(1) the condition for which nontrivial solutions for the 
value exist in the system 

(2) the condition that these o r l values are such that the 
motion be stable, i.e., that the real parts of all the CT t values 
be less than zero. 

In order to fulfill condition (1) above, we substitute 
expressions (Lid) into equations (ii3) through (b7) and realizing 
that the principle of superposition of solutions holds for 
linear differential equations, drop the summation notation and 
obtain the following relations: 

1. Barge 



a) Lateral Force 




049) 



b) Yawing Moment 




( 50 ) 



2. Towboat 



a) Lateral Force 



- 15 




b) Yawing Moment 
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c) Rudder 



\ = % ^ 
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(53) 



We now have a system of linear homogeneous equations in 
>/' t , , t , £ , and S t . The condition that a nontrivial solu- 

tion exists among them is that the determinant of their coef- 
ficients be identically equal to zero. This means that 
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Expan. ion of this determinant leads to the following rela- 
tion: 

* «-j 

1 C . cr - 0 . 

j-o J (60) 

wh r : 

C Q = (ah - /d)(AH - LD) 

C 1 = (ah - id) (AJ + BH - LF - MD) + (aj + bh -/f - md)(AH - LD) 

C 2 = (ah - /d)(BJ ♦ EH ♦ LG - MF - ND) 

+ (aj + bh -if - md ) ( A J + BH - LF - MD) 

4 (bj + eh - ig - mf - jg)(AH - LD) 

C 3 = (ah - l d)(EJ + Mg - NF) 

+ (aj + bh - if - md)(BJ + EH + Lg - MF - NT)) 

4 (bj + eh -ig - mf - jg) (AJ + BH - LF - MD) (6l) 

4 (ej - mg - jf ) (AH - LD) 

-- (.ah - id) Kg + (aj + bh - £.£ - md)(EJ + Mg - NF) 

+ (bj 4 eh - l g - mf - jg)(BJ + EH + Lg - MF - ND) 

4 (ej - mg - jf)(AJ + BH - LF - MD) 

- jg (AH - LD) - g 1 1 L 
= (aj + bh ~if - md) Kg 
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To fulfill condition (2) on page (lh), we need look to 
the theory of equations as developed for polynomials. In Ref. 
3, we find that the conditions for stability (the real parts 
of all the roots are less than zero) are stipulated by the 
following considerations: 



(1) all C . >0 
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(2) all deternintal relations obtained from 
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first the first two rows and columns 
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second the first three rows and columns 



l. e. 



J 0 



0 

C 1 

C, 



(6M 



> 0 



( 65 ) 



( 66 ) 



(67) 



ir* i o or until ill u^h « Apre.vir r.r ha/< b* ^ r* wr M n 1 r. - 
iUding (6 f j . 

Tn< . < inojualiti- , when aatirfiod, indicate ^tabi iiy 
r or the* -y tern unier consideration. For any t wo \ articular 
bodi* , a towline length and towpoint location^ on both the 
barge and towboat may be ietermined to give desirable tab! ity 
haracteri sties. This, however, requires considerable nur.eri- 
ai anaiy i. a~ well as time, unless suitable co piling fa- 
ilities are available. 

The development of the theory as presented here io con- 
idered to be a basis for future numerical investigation of 
the directional stability of any two bodies separated by a tow- 
line of arbitrary length attached to arbitrary points along 
the centerline of each body. 



A SIMPLIFIED 1 Hfc’ORY 



The foregoing nrvaly i ' ha- boon dfrvotel to the dir' ction- 
al stability of a tow -d configuration of two bodi-_, o-parated 
by a towline of arbitrary length, attached at arbitrary points 
along the centerline of each of the bodies, and further con- 
fined by a control rule for the towing body proportional to 
its heading deviation. 

If is evident from this analysis that the results ob- 
tained defy, without considerable numerical calculation, ra- 
tional qualitative results. Thus in order to correlate quali- 
tative as well as some quantitative predictions with data from 
test measurements, certain rather broad simplifications will 
be made. 

It will be assumed that 

(1) the towline is horizontal, and taut (as before): 

(2) the barge does not roll or pitch 

(3) the towboat proceeds on a straight and prescribed 
course regardless of the motion of the barge. 

With this latter assumption, all the towboat functions 
vanish and only the barge functions remain. 

It is also convenient to define the transformation and 
dimensionless coefficients presented on pages 9 and 10 as: 
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With these in mind, equations (h3) and (hh) reduce to: 
Lateral Force 
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Yawing Moment 
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Reasoning similar to that presented on pages 13 through 17 
leads to the stability equation 
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The motion will be stable if and only if all the Cj 
values are positive and the following relations are satisfied: 



- Cl C 2 - C 0 c 3 > 0 



( 73 ) 



and 






0 0 



* (C 1 C 2V C 1 2 V C 0 C 3 2)>0 



( 71 .) 



In the case of the untowed body, it is noted (Ref. 1) 



that 



/ / — — 
-n m 2 CTj ct 2 



( 75 ) 



N'i Y'p + N' fi 

n'y'p - m 2 N 'jj = - n ' m 2 (*1 + 

where & { and & 2 are the dynamic stability indices of the un- 
towed body. Then the C . values may be written as 

r _ ' / 

L 0 ~ n m 2 

C 1 = ~ n> *2 (°l + ° 2 ) 

C 2 ' f'p T *! *-p-> ' n 'p* "2 5 | 5 I_ 

c 3 ■ [«; ''V > t'* o'i— ; * 

c i. ■ [V 

Without further mathematical manipulation, the following 



( 76 ) 
T 

f 



analysis may be made: by definition, all the values must 

be positive. Then: 



1) Cq > 0 always 



- 21 



2) C x > 0 for: 

a) all dynamically stable untowed bodies, since 

a j and a < 0 ; 

b) all dynamically neutrally stable untowed bodies, 
since 

a x = 0 and ct 2 < 0 ; 

c) all dynamically unstable untowed bodies for 
which 

|<7,| < |cr 2 | ( where o r l > 0 and & 2 < 0 . 

As a consequences of this, a towed body cannot 
be stable on course if I® 1 1 1 > I c’" 2 1 and 

l^i I > 0. 

3) C 2 > 0 for: 

a) all dynamically stable untowed bodies, since 
the towing functions are positive and since 
both CTj and o 2 are negative; 

b) all dynamically neutrally stable untowed bodies, 
since 



o j = 0 and & 2 < 0 ; 

c) all dynamically unstable bodies for which 
x p T' n 2 (1 + £ ) + n + n ' m/ o, a 2 > 0 . 



M 



> 0 for all towed bodies 
since by definition, (* p / Y‘ 



ally (y'+ 

r > 0 . 



*2 ) > 0 * 



(see item 5 below), 
p — N p) > 0, and usu- 
S ' ^ >0 and T / and 



1 



' 0 . 



) C > 0 by definition. Jincc 



* * ' » - v 



> b 

* c- 



— \ 



this condition may be fulfilled either by increa ing 
x ' (moving the towpoint farther forward), or, if 
this is in a practical sense impossible, by increas- 
ing the drag in such a way that the static moment 
rate will remain unchanged or is decreased. The 
condition Cj > 0 defines the location of the tow- 
point: it must always be a distance x'p forward of 

the center of gravity of the ship, which is greater 
than the distance between the center of gravity and 
the center of pressure of the static force Y pG. 

Assuming that the C. values are all positive, it is ne- 

J 

cessary to examine the determinantal relationships 

a) A, =- n'n 2 ' (5 x + ) £x p ' T ' (1 + JL ) m' + n ' y, + n ^2 ] 

-"X {[*p (l +X f)T'+ [*p ) 

*\v 



n'm' (a, a, * 2 | x' m/ (1 ^ 

’l-m' + *„' )-\ {J-Lj> 



' v 1 \ ' 

x r ] p-* 0 



1 r / - 


+ 


- ?d 


- 



X P + m 2 *| 



-n'm/ (cr, +cr 2 ) cr, o 2 - t' (a x *cr 2 ) Q + ir > 0 



where 



y = X ' rn ' (l + X JP ) + JL 

v /) 2 7 2' 



and 




and Q, It, and n ' m' 2 must be positive. The following cases 
will now be exanined: 

i) Untowed body, dynamically stable. Here, 

- n m 2 (a l + B 2 ) (7 , <7 2 > 0 

always for 5", ~a 2 > 0 and (& x * & 2 ) < 0. Thus, the 
sense of the inequality depends on the magnitude of 
t'H’, hence, to achieve stability, it is desirable to 
keep this term small. This may be accomplished by 
proper selection of , the towpoint, and „ the 
towline length. 

ii) Untowed body, neutrally stable. The first term 
vanishes since a 2 =0; thus H must be made small 
in order to satisfy the inequality, and insure sta- 
bility. To accomplish this, the same reasoning is 
applied as in a). 

iii) Untowed body, dynamically unstable. Here, ^ > 0 
and as long as o 2 < 0 and \a l \ < |a 2 then 

-n m 2 (a l + a 2 ) a, a 2 < 0 

Thus, for stability, it is desirable to make fl as 
small as possible and to increase the drag in such a 
fashion that the inequality is satisfied. 






I 
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b) a 2 - Cl c 2 c 3 - c x 2 c u - c Q c 3 z > 0 



^3 (C 1 C 2 * W 



'1 "h 



C, > 0 



c 2 c. 

a 1 h 

A > CT“ 



> o 



(76) 



If it is assumed that the conditions specified for A, 
have been met , then the remaining factor to consider is the 
magnitude of C^ 2 (C^/C 3 ). C-^ is determined by the values of 

a 2 and n' i* 2 , and 0 << C^/C^ << 1. It is therefore de- 
sirable to make as small as possible, still keeping it 
positive, and thus insure the stability of the system. 

The criteria of directional stability will now be deter- 
mined for the case where the towline has increased in length 
to a point where the terms containing its reciprocal are 
negligible. Mathematically, this is achieved by letting the 
towline length increase without limit. In this case, it is 
found that 



A 2 



/ / 

- n m 9 




+ 



fo-j + 7 2 )(x' p T'+ no j cr 2 ) 

*'p *' / *-* / / 3 ] 7 ' / ] >0 



Since n m 2 , T / and (x p Y - K ' e ) are positive, this ex- 
pression will be positive if 

(<"/ (o, ♦ o 2 )(x' r'+ no, o 2 ) ♦ (x' p Y' fi - v \) t' <] o 



or 
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(x 



/ 

p 



/ 




V 




/ 

< 



m 



/ 



2 



/ / _ 
T + n cr x 



o 2 ). 



(* l + *2* (*p 






I 

I 



For all dynamically stable ships, the right-hand side of the in 
equality is positive; thus (x f) ' Y ^ — N / p)T ' must be kept as 
small as possible if the resulting motion is to be direction- 
ally stable. For dynamically unstable ships, the directional 
stability depends on (x ' T +n' a, > 0. Since in this 

case a 2 <0, it must be specified that for directional sta- 

/ ' i — — 

bility, x p T> — na , l ar J i n addition to the condition that 
( x p V'p -N V T'< -m' (5 X ♦ a 2 ) (x' 7 n'a, * 2 ) 

To investigate the consequences of decreasing the towline 

length to zero, i.e., either towing the ship fron its bow or 

with a bridle, we discard the terms containing the towline 

length 2 in equations ( 69 ) and (70) on page 22 since they have 

no significance in this situation. The resulting stability 

equation is then a cubic of the form of 
3-j __ 

2 C . cr = 0 (80) 

j = 0 J 1 

where 



n m 0 








("V* - 


/ ,, / . ) = / f 
m 2 /V^/ — n m 2 


(cr x 


cs 

lb 


. / t 


/ , i / / / 


/ 


/ / ~ 


< Y j, -m l 


) N 0 ~ Y ^ * +m 2 X p 


7 = 


n m 2 cr 



C3 - (*p Y' fi -N' p )T' (61) 



The criteria for stability specifies that all the C. 

J 

values must be positive and that the following determinant is 



satisfied : 












I 



I 
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C 



1 




A c 



3 



( 02 ) 



0 



0 



C 



3 



This reduces to (C^ - Cq C^) >0 or 

) | n a, a 2 > x' p 7 'J— n'm/ (x' p y ' 0 V ^) 7 '> 0 




Hence, the motion will be stable if 




/ 



( 83 ) 



For dynamically stable ships, the left-hand side of the inequal- 
ity is positive; consequently, the stability or instability of 
the towed ships is essentially a matter of towpoint location. 

For dynamically unstable ships, however, the extent of 
instability will be an important factor in determining the di- 
rectional stability, since the sign of the left-hand side of 
the inequality will depend on the sign of (n'a l cr 2 + x p T )■ 
Thus, for the resulting motion to be stable, ("V, cr 2 + x p ^ ^0. 
Hence, the dynamic stability criterion must be small if this 
condition is to be satisfied. 

Reducing towline length, however, is of very little prac- 
tical consequence, since an excessive decrease in towline length 
causes hydrodynamic interaction between the towing and towed 
ships, a fact that will cause variation in the hydrodynamic 
"constants". Thus, the basic assumptions of the analysis break 



down and the conclusions derived become untrustworthy. 
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KERICAL ANALYSIS 



The results of an investigation to detcrnine the hydrody- 
namic coefficients of a typical barge form with various skeg 
configurations made on the rotating arm of the Experimental 
Towing Tank afford the data necessary for an analysis based on 
the foregoing treatment. 

The particular model considered had been found to be 
stable according to the free-towing standards set by the Experi- 
mental Towing Tank. These standards stipulate that if a barge 
is towed from its bow with a towline length of twice the barge 
length, that the barge should not be laterally displaced more 

■f ^ 

than -5 of its beam. The towing device is a carriage that pro- 
ceeds along the tank rail at constant speed for the major portion 
of its travel. 

The accompanying sketch shows the bare hull of the barge 
model together with the skeg modifications investigated on the 
rotating arm (Ref. 6). The accompanying tabulation gives de- 
tails of the configurations as well as the comparable results 
of the standard free-towing stability criterion set by the Ex- 
perimental Towing Tank where these results are applicable. 




Model Particulars: 
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n' - .00219 
l = 6.5 ft. 
p = 1.9b 



0.0351 



^ = o.5oa 



", - 0.0219 
?' = 2.0 



Case I - bare hull (unsatisfactory by free-towing method) 

Case II - bare hull and skegs (5° to £ ) (unsatisfactory by 
free-towing method) 

Case III - bare hull, skegs and stabilizing prisms (satisfac- 
tory by free-towing method) 

Case IV - bare hull and slotted skegs (not tested by free- 
towing method) 

Case V - bare hull and slotted skegs deflected to 6° (not 
tested by free-towing method) 

Case VI - bare hull and slotted skegs deflected to 12° (not 
tested by free-towing method) 

In each case the model was tested on the rotating arm at 

y 

a speed length ratio of = O.I 163 where V is measured in 

knots. The following values of the hydrodynamic coefficients 



were determined: 



Case 


v' 

Y P 


j v / 

A 0 


l / = t' 


v ' - 

r <J> 




I 


0.008 


0.010 


- 0.00250 


0.0038 


-0.0020 


II 


0.020 


0.0057 


-0.00280 


0.00b9 


- 0.0029 


III 


0.020 


0.00b9 


-0.0032b 


0.0055 


- 0 . 003 b 


IV 


0.022 


0.0069 


-0.00280 


0 0038 


-0.0028 


V 


0.020 


0.0060 


-0.00280 


0.0037 


-0.0029 


VI 


0.019 


0.0058 


-0.00280 


0.0029 


-0.0030 



Calculation in accordance with the modified development 
(pg.21) and Ref. 1 leads to the following values for the untowed 
stability criteria ( a t ) togeth- r with the towed stability 



criteria a. w tl a and \ . 

0 





CASE I 


CASE LI 


CASE IJJ 


case i_y 


case v_ 


CASE Y1 


a l 


♦ 1 047 


>0 237 


-*0 0735 


*0 363 


♦0 303 


♦ 0 3 1 1 




2. i** 


2 131 


-2 1*5 


-2 26* 


-2 197 


-2 222 


Co X 10 4 


0 76*7 


0 76*7 


0 7 ’7 


0.76*7 


0 76*7 


0 76*7 


C, x I0 4 


0*772 


l 4559 


1 6314 


1 4646 ; 


l 455* 


1 45 Q 1 


C 2 x 10 4 


-1 06 


1 270 


0.630 


0 0239 


0 1448 


0 1248 


C 3 x I0 4 


-0.0532 


0 32 5 Q 


0.4227 


0.31 12 


0 30*2 


0 2*34 


O 

X 

O 


-0 0742 


0 0624 


0.0*52 


0.0602 


0 05*8 


0 0546 


‘ l x 10 H 




0. 1426 


0.702* 


-0.2042 


-0.0261 


-0 0422 


A 2 x lO 14 




- 0 . 095 * 


0.0704 


-0. 1926 


-0. 1326 


-0 1302 




0.82 


~ 0.576 


[-0.031 1 


[ 0.1254 


0.0*05 


" 0 0919 


a 2 


r 0.044 


_± 0. 4775 i 


I to. 509 i 


[± 0 . 4606 i 


± 0 . 4634 i 


_ ±0.45071 




|± 0 . 248 1 


-0. 1^35 


-0.2405 


-0 1735 


-0. 1800 


-0. 1744 




-1 873 


-1 8155 


- 1 8195 


- 1.9825 


-1.895 


-1 .9200 



The tabulation shows that none of the cases considered is 



dynamically stable since in each case one of the stability 
roots (5^) is positive. Of all the skeg configurations, only 
Case III produced a directionally stable barge when towed 
(since this is the only condition where all of the real parts 
of the values are negative). Case III was also considered 
stable as a result of the free-tow T ing tests previously des- 
cribed. However, the drag of the barge in this case is notice- 
ably higher than for any of the others. 

It will be seen that none of the other skeg configura- 
tions creates such high resistance. While these other cases 
are not directionally stable, it appears that the instability 
may well be reduced without creating excessive drag. It does 
not seem too unreasonable to assume that further experimental 



investigation will yield .;kcg foms that will produce direction 
ally stable barges without adverse drag increments. 

The f resences of complex values (paired in the table) 
indicates that the motion is characterized by oscillatory modes 
In all cases except Case III, the oscillation is undamped and 
tends to build up with time. This fact imposes continuous al- 
ternate yawing on the towed ship. In Case III, which is di- 
rectionally stable, the oscillations still appear but are 
damped and tend to diminish with time. 



CONCLUSIONS 



The preceding analysis shows that while the stability cri- 
teria and equation, for the case when two bodies are considered, 
have been determined, their evaluation for a given set of physi- 
cal and hydrodynamic constants is extremely laborious. It is 
difficult, if not impossible, to gain an extensive amount of 
qualitative information from the coefficients of the stability 
equation and the stability criteria because of the complicated 
way in which the various terms comprising these expressions ap- 
pear. 

In the simplified analysis where only the towed body is 
considered the following generalizations may be 'made. 

1. If a towed ship is to be directionally stable, the tow- 
point must be located forward of its center of gravity greater 
than the distance between the center of gravity and the center 

of pressure of the static lateral force. 

2. Depending upon its degree of stability, mass, towpoint, 
resistance, and moment of inertia, a dynamically stable ship may 
become unstable when towed with a long enough towline. 

3- Depending upon its towpoint, mass, moment of inertia, 
and damping functions, a dynamically unstable ship may become 
directionally stable when towed with a short enough towline, pro- 
viding that there is negligible hydrodynamic interaction between 
it and the towing ship. 

b. The motion of a towed body, which may be either di- 



rectionally table or unstable, i oft^n characterized by o di- 
latory mode^. 

5 . In cases where a ship is extremely instable dynami- 
cally, no possible towing configuration can make it direction- 
ally stable . 

6. Increase in hull drag aft causes variation in each of 
the hydrodynamic parameters. Sufficient drag will cause most 
towed ships to become directionally stable. The experimental 
data presented indicate that the arag increments nay be mini- 
mized and directional stability achieved by alteration and im- 
provement of skeg conf igurations . 



A FPL’.' X 
IJOMLNt LATURE 



All lin ar dinr n ion are in fo> t ani angjlar dimensions 
in radians. The tine unit is th recond. 
t barge length 

x p barge towline point of attachment forward of barge 
center of gravity 

CG barge center of gravity 

towboat length 

* p towboat towline point of attachment aft of towboat 
center of gravity 

CG towboat center of gravity 

0 barge heading angle 

$ barge yaw angle 

0 towboat heading angle 

towboat yaw angle 

S towboat rudder angle 

— • • 

S towboat rudder angle order 
t time 

7 towboat rudder time lag 

V barge linear velocity 

V towboat linear velocity 

2 tow line length 

L over-all length of tow, i.e., from bow of towboat to 
stern of barge 

V mass of barge 

V mass of towboat 

V, longitudinal virtual na^s of barge 
M lateral virtual mass of barge 
M, longitudinal virtual ma s of towboat 
\l 2 lateral virtual mass of towboat 
/ moment of inertia of barge 

1 B monent of inertia of towboat 

/ virtual ncm- nt of inertia of ba^ge 



/ virtual mom»nt of inertia of towboat 
l w moment of inertia of displace l fluid of bargo 

mom- nt of inertia of dirplaced fluid of towboat 

k , longitudinal virtual ma. s coefficient of barge 
fc 2 lateral virtual mass coefficient of barge 
k j longitudinal virtual mass coefficient of towboat 
k 2 lateral virtual mass coefficient of towboat 
k' virtual inertia coefficient of barge 

T' virtual inertia coefficient of towboat 

/ towline tension 

p density of water 

} hydrodynamic lateral force of barge 

/V hydrodynamic lateral moment of barge 

Y hydrodynamic lateral force of tow'boat 

N hydrodynamic lateral moment of towboat 

first, second, and third derivatives with respect 

to time 

/ indication of dimensionless coefficient 

subscript notation denotes a partial derivative 
with respect to the subscript; thus 

ik = Y e ^ = ^ etc ' 

m j dimensionless longitudinal virtual mass coefficient 
of barge 

dimensionless lateral virtual mass coefficient of 
barge 

m x dimensionless longitudinal virtual mass coefficient 
of towboat 

m 2 dimensionless lateral virtual ma s coefficient of 
towboat 

n' dimensionless virtual moment of inertia coefficient 
of barge 

n' dimensionless virtual moment of inertia coefficient 
of towbo* 

All o*her auxiliary cynbols are defined where they appear 
in the text. The exact expressions for the dimensionless factors 
are given in each instance where they are first introduced. 
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